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The electronic structure and magnetic properties of CaCrOs have been calculated by two methods, 
including hybrid-exchange density-function theory and density-functional theory + U. The com¬ 
puted densities of states from both of these methods are in a qualitative agreement with the previous 
x-ray spectroscopy. On the other hand, the opening of the band gap separates them apart, hybrid- 
exchange density-functional theory always gives a finite band gap, down to ~ 1.2 eV from HSE06 
functional, whereas by tuning the Hubbard-17 parameter down to 0.5 eV, a conducting state with 
AFM-C (defined in the text) spin configuration can be achieved. From hybrid density-functional 
theory, the computed nearest-neighbouring exchange interaction along the c-axis and in the a6-plane 
are ~ 4 meV and ~ 6 meV (anti-ferromagnetic), respectively, which are qualitatively in agreement 
with the previous magnetic measurements. The density-functional theory + U has predicted a sim¬ 
ilar series of exchange interactions. These anti-ferromagnetic exchange interaction, together with 
the in-plane anti-ferro-orbital ordering will induce a spin-orbital frustration, which could play a role 
for the abnormal electronic properties in CaCrOs. In hybrid-exchange density-functional theory, an 
abrupt reduction (~ 0.2 eV) of the majority-spin band gap of the ferromagnetic state between 60 
K and 100 K has been found as lowering temperature, which shows a strong link to the previous 
optical conductivity measurements in [A. C. Komarek, et. al., Phys. Rev. B 84, 125114 (2011)]. 

In sharp contrast, the density-functional theory + U methods predicted AFM-C state as the lowest 
AFM state for the crystal structure measured below 90 K, above which AFM-A is however the 
lowest. The closely related concepts including electron-hole liquid and surface-plasmon-mediating 
spin-spin interactions have been discussed as well. 

PACS numbers: 75.47.Lx, 71.15.Mb, 75.30.Et, 75.25.Dk, 75.25.-j, 73.20.Mf 


I. INTRODUCTION 


Oxides, widely known as ceramics in the modern world, 
can be even found in the ancient pottery. Among these, 
transition-metals oxides (TMOs) have copious interest¬ 
ing physical phenomena, ranging from high-transition- 
temperature superconductors 11], colossal magnetoresis¬ 
tance [2|, and multiferroicity |3] to plasmonics Q. For 
example, chromium dioxide (Cr 02 ), which is a half metal 
(conducting in one spin channel but insulating in the 
other), has a great potential for spintronics @, thus ex¬ 
ploring electron spins to store and transport information. 
These fascinating physical phenomena, as well as their 
important applications in the development of new tech¬ 
nologies, hugely depend on the understanding of the in¬ 
teraction between (i) the localised d-electron of transi¬ 
tion metals and p-electron of oxygen and (ii) between 
the spin and orbital degrees of freedom. Moreover, cova¬ 
lent, metallic, and ionic bonds can all be found in TMOs, 
hence leading to possible boundary crossing between dif¬ 
ferent electronic and magnetic states. The underlying 
physics behind these interactions are also manifested in 
the phenomena of fundamental interest, such as spin- and 
orbital-orderings. 

Cr 4+ -based perovskites, such as CaCrOs, SrCrC> 3 , and 
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FIG. 1: (Colour on line.) The conventional cell of perovskite 
CaCrOs is shown. Ca is depicted as green ball, O as red, and 
Cr as yellow. 


PbCrOs, have attracted considerable attentions recently 
owing to their peculiar electronic properties. Es¬ 
pecially in CaCrOs (whose crystal structure is shown 
in FiglU), the complicated interplay between magnetism, 
conductivities, and lattice vibrations has been the major 
obstacle for a clear observation of its electronic ground 
state E0- The question, whether CaCrOs is metal¬ 
lic or insulating, is still under active debate. A com- 
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prehensive experimental study has suggested that the 
anomalous properties of CaCr0 3 are related to the Cr-0 
bond instabilities owing to the spin and orbital orderings, 
which could result in a crossover from localised to itin¬ 
erant behaviour j6j. However, a more resent study, com¬ 
bining experiments and first-principles calculations, has 
shown that CaCr0 3 is a rare case of anti-ferromagnetic 
metal owing to the co-existence of anti-ferromagnetism 
and metallic-like conductivity 0 , The theoretical 
work therein replying on local spin density approxi¬ 
mation (LSDA) + U supported that the c-type anti¬ 
ferromagnetic (AFM) state, in which the spins are in fer¬ 
romagnetic (FM) ordering along the c-axis but AFM in 
the a6-plane, is the ground state @]. the electronic state 
was predicted to be metallic within LSDA, but insulating 
once an appropriate value of XJ was included. Another set 
of first-principles calculations 0 based on generalised 
gradient approximation (GGA) + U have suggested that 
the c-type AFM is a metallic state with an intermediate 
U. To the author’s knowledge, most of the previous the¬ 
oretical work has been done either within pure density- 
functional theory (DFT) or DFT + U [§, 0 ] for only one 
crystal structure. Against this background, Komarek, et. 
ah, have measured the crystal structures in a wide range 
of temperatures from T = 3.5K to 300A'. This work has 
provided a good opportunity to study the behavior of the 
electronic structures as a function of temperature, which 
is yet to be studied thoroughly. 


In this paper, the electronic structure of CaCr0 3 was 
computed within two theoretical frameworks: (i) hybrid- 
exchange density-functional theory (HDFT), PBEO 0 
and HSE06 [l2j], both of which are hybrid-exchange den¬ 
sity functional free of any adjustable parameter, and (ii) 
DFT + U methods where the Hubbard-f/ can be ad¬ 
justed, according to different crystal structure reported 
in ii- The total energies of ferromagnetic (FM), 
A-type AFM (AFM-A, AFM along the c-axis but FM 
in the ab plane), C-type AFM (AFM-C), and G-type 
AFM (Neel state, AFM-G) states were computed care¬ 
fully. The nearest-neighbouring (NN) exchange interac¬ 
tions along c-axis and in the afo-plane were then extracted 
by comparing the total energies of different magnetic 
states. We have found that AFM-G spin configuration 
is the ground state, which is different from that observed 
and predicted in Ref. Q ■ However, the total energy dif¬ 
ference between the AFM-G and AFM-C states is < 30 
meV within the reach of room-temperature. These calcu¬ 
lations suggest that spin- and orbital- orderings co-exist 
in CaCr0 3 . The rest of the discussion is organised as the 
following: in m the computational details are given, in 
m the calculation results are presented and discussed, 
and in HIVI some more general conclusions are drawn. 


II. COMPUTATIONAL DETAILS 

A. Hybrid-exchange density-functional theory 

Calculations for the electronic structures and magnetic 
properties of Cr02 and CaCr0 3 were carried out by using 
DFT and a variety of approximate exchange-correlation 
functionals, including PBEO 0 , HSE06 01. and HISS 
0 ], implemented in the CRYSTAL14 codefla]. The ex¬ 
perimentally determined latti ce p arameters and atomic 
coordinates are adopted here 0 , 0 . The basis sets of 
Ca 0 , Cr 0 , and O 0 ] that are developed for solid- 
state calculations were used. The Monkhorst-Pack sam¬ 
plings [20] of reciprocal space are carried out choosing a 
grid of shrinking factor to be 7 x 7 x 5 to be consistent 
with the ratios among reciprocal lattice parameters. The 
truncation of the Coulomb and exchange series in direct 
space is controlled by setting the Gaussian overlap tol¬ 
erance criteria to 10~ 6 ,10 -6 ,10 -6 ,10 -6 , and 10 -12 0 ]. 
The self-consistent field (SCF) procedure is converged to 
a tolerance of 10 -6 a.u. per unit cell (p.u.c). To ac¬ 
celerate convergence of the SCF process, all calculations 
have been performed adopting a linear mixing of Fock 
matrices by 30%. 

Electronic exchange and correlation have been 
described using the PBEO hybrid-exchange, sliort- 
(HSE06), and middle (HlSS)-range-corrected function¬ 
als. Apart from a partial elimination of the self¬ 
interaction error, these functionals balance the tenden¬ 
cies to delocalize and localize wave-functions by mixing 
a quarter of Fock exchange with that from a generalized 
gradient approximation (GGA) exchange functional 0 - 
The performance of the hybrid functional, e.g., B3LYP or 
PBEO has previously been shown to provide an accurate 
description of the electronic structure and magnetic prop- 
erties for both inorganic and organic compounds [221 . |23I | . 

B. Density-functional theory + U 

The DFT + U method developed in the Quantum 
Espresso code [0 has also been used for a benchmark 
and further comparison with HDFT method. The PBE 
(Perdew-Burke-Ernzerhof) exchange-correlation func¬ 
tional 0, along with a range of values of Hubbard- 
U for Cr, has been chosen. In this method, the same 
set of lattice parameters and atomic coordinates as in 
HDFT calculations are used for consistency. The Van¬ 
derbilt ultrasoft pseudoptentials developed for the PBE 
exchange-correlation functional have been used for all the 
elements throughout DFT + U calculations. The same 
Monkhorst-Pack sampling has been used for DFT + U as 
the HDFT calculations. The planeO-wave cutoff energies 
of ~ 640 eV and the threshold of SCF energy convergence 
of 10~ 6 eV were employed. A linear mix of 30 % of the 
Fock matrix has been used to accelerate the SCF conver¬ 
gence. The electron correlation associated with 3d states 
of Cr are described by DFT + U methods. A range of val- 
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ues of U = 0.5,1.0, 2.0 and 4.0 have been chosen to show 
the transition between metallic and insulating states. 

C. spin Hamiltonian 

The Heisenberg model [2l| for the spin-spin interac¬ 
tions in CaCrC >3 is defined here as, 

H = J c ^2 Si ■ Sj + Jab ^2 Si ■ Sj , ( 1 ) 

ijdzC ij£ab 

where J c is the exchange interaction between the NN 
spins along the c-axis and J a t between the NN spins in 
the ad plane respectively. The i and j label the Cr sites. 
The exchange interactions are determined by 

Jc = {E FM ~E a )/8S 2 , ( 2 ) 

Jab = {E FM ~E C )/8S 2 , (3) 

where E F m, E a , E c , and E g are the total energies for a 
conventional cell in the FM, AFM-A, AFM-C, and AFM- 
G states, respectively. S' = 1 is adopted for a Cr 4+ ion. 

III. RESULTS AND DISCUSSIONS 
A. HDFT calculations 

1. Electronic structure of CaCiE)s 

The projected density of states (PDOS) onto the 3 dr 
orbitals of Cr (the first row), 2p-orbitals of O (the second 
row), and spin densities (the third row) of the AFM-A, 
AFM-C, and AFM-G states, calculated by using PBEO 
for the crystal structure reported in Ref. [131] . are shown 
in Fig]2] The electronic states of FM, AFM-A, AFM-C, 
and AFM-G were predicted to be insulating with a band 
gap ~ 2 eV. The band gaps in different AFM spin config¬ 
urations, computed by using three functionals, are listed 
in Tabled] The lowest band gap is given by HSE06 in the 
AFM-C state, which is ~ 1.44 eV for the crystal struc¬ 
ture reported in Ref. fl3|, while for the crystal structure 
in Ref. [Hi, the smallest band gap is ~ 1.2 eV. Such nar¬ 
row semiconducting band gap could result in high con¬ 
ductivities and a facile metal-insulator transition. How¬ 
ever, PBEO and HISS have shown slightly larger band 
gaps, up to ~ 3 eV. The AFM-G state is predicted to 
be the lowest state; this is different from the previous 
theoretical and experimental work [jj. This is probably 
owing to the larger Hubbard-Z7 embedded in the hybrid- 
exchange-functional methods, as compared to the value 
of U adopted in the previous local-density approxima¬ 
tion + U method 0- However, the difference between 
these three AFM states is up to ~ 38 meV, which is 
close to the energy scale of room temperature. Espe¬ 
cially the smallest computed energy difference between 
the AFM-A and AFM-C states is ~ 2 meV. Therefore, 


the stability of these magnetic states is fragile respected 
to thermal fluctuation and lattice vibration. Thermal 
and/or lattice fluctuations will further induce a magnetic 
disorder, which could prevent a proper observation of the 
electronic state in CaCrC> 3 . The 3d-orbital PDOS near 
VBM is dominated by that from d x 2 _ y 2 , d xz , and d yz . 
The 2p-orbital PDOS of O in the AFM spin configura¬ 
tions share the similar feature; one significant peak is 
near the VBM with the other close to —3.7 eV. The com¬ 
parison between the PDOS of 3d-orbital and 2p-orbital 
suggests there could exist a strong hybridisation between 
them, which is especially strong near the VBM. This is 
in agreement with the claim made in the Ref. 3] that the 
electronic properties of CaCr 03 is driven by the strong 
2p-3d hybridisation. This strong bonding between Cr 
3d-orbitals and O 2p-orbitals might be able to shed some 
light on the high conductivities of CaCr 03 . Moreover, 
the O 2 p PDOS is more dominant than the Cr 3d PDOS, 
which is qualitatively different from Cr02- This could 
suggest that the conductivity in CaCr 03 originates from 
predominantly 2p-orbitals rather than 3d. 

As shown in the spin density of the AFM-A state 
(Figl2k), the spins are anti-aligned along the c-axis, but 
aligned in the ad-plane. In AFM-C (Fig[5Jr), the spins are 
aligned along the c-axis, but anti-aligned in the ad-plane. 
And in AFM-G, a neel state is present, i.e., anti-aligned 
both along the c-axis and in the ad-plane. Although there 
is a clear difference between these magnetic states, they 
share a common feature that there exists an anti-ferro- 
orbital ordering (illustrated by yellow arrows) in the ad- 
plane, whereas a ferro-orbital ordering is formed along 
the c-axis. This will result in a frustration between spin- 
and orbital-orderings in the ad-plane for all the AFM 
states. People have shown in in iron-based superconduc¬ 
tors that spin-orbital frustrations that are intrinsically 
due to quantum fluctuations can induce abnormal metal 
state [26j . 

We have also analysed the spin densities projected to d- 
orbitals. As shown in Table ID the spin densities are dom¬ 
inated by d xz , d yz , and d x 2_ y 2 orbitals (~ 1.7 ps for these 
three orbitals as expected for spin-1) from all the func¬ 
tionals used here. The d x 2_ y 2 is occupied mainly because 
the tetragonal distortion, which will lower d x 2_ y 2 and lift 
up the d z 2 orbitals. The PDOS calculated here (Fig]2|) 
is in a good agreement with this physical picture. The 
nearly even mixture of d xz and d yz is in agreement with 
the previous theoretical results that claimed a ^-active 
orbital ordering Q. The two d-orbitals, d xz and d yz , are 
responsible for the complex orbital ordering in the ad- 
plane as shown in Fig[5^, h, and i. The orbital ordering 
in the ad-plane is formed by the alternating patterns of 
the ^=(| d xz ) + | d yz )) and -E(\d xz ) - | d yz )) orbitals. 

Komarek, et. al. 0 have observed the evolution of 
the crystal structure of CaCr 03 when changing tempera¬ 
tures. The corresponding electronic structures have also 
been computed by using HSE06 hybrid functional. The 
DFT band gaps of the majority spins in the FM config¬ 
uration as a function of temperatures have been plotted 
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FIG. 2: (Colour on line.) The PDOS of Cr d-orbitals, O 2 p orbitals, and spin densities on Cr, computed by using PBEO for 
the crystal structure reported in Ref. El, are shown for A-type, C-type, and G-type AFM spin configurations for CaCrOs in 
columns. The PDOS for d z 2 is depicted in red, d xz in green, d yz in blue, d x 2_ y 2 in cyan, and d xy in black. The PDOS for p x 
is depicted in red, p y in blue, and p z in black. Zero-energy is chosen to be at the valence band maximum. In spin densities, 
spin-up is in red and spin-down in green. Notice that the orbitals are anti-ferro-ordered in the a6-plane (the orientation of 
orbital is illustrated by yellow arrows). 
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PBEO HSE06 HISS 

AFM-A AFM-C AFM-G AFM-A AFM-C AFM-G AFM-A AFM-C AFM-G 


EgapieV ) 

2.65 

2.15 

2.50 

1.96 

1.44 

1.89 

2.96 

2.49 

2.97 

d z 2 

0.11 

0.11 

0.11 

0.11 

0.11 

0.11 

0.10 

0.10 

0.10 

dxz 

0.52 

0.55 

0.51 

0.51 

0.55 

0.52 

0.50 

0.52 

0.49 

dy Z 

0.54 

0.53 

0.56 

0.55 

0.53 

0.54 

0.54 

0.52 

0.54 

dx 2 —y 2 

0.60 

0.69 

0.51 

0.70 

0.69 

0.70 

0.72 

0.71 

0.72 

dxy 

0.13 

0.14 

0.14 

0.14 

0.14 

0.14 

0.13 

0.13 

0.13 


TABLE I: The projected spin densities onto Cr d-orbitals, computed by using PBEO, HSE06, and HISS are shown for the 
AFM-A, AFM-C, and AFM-G states. 


in FiglHI The band gap is indirect; the valence band 
maximum is at the T point whereas the conduction band 
minimum (1/2,1/2,0). Notice that there is an abrupt 
jump of ~ 0.1 eV in the FM band gap between 90 K 
and 100 K, which coincides with the jump in the mea¬ 
surements of resistivity and optical conductivity in Ref. 
0. In contrast, the minority band gap in FM configu¬ 
ration maintains as ~ 2.7 eV. The other band gaps have 
not been affected by the temperatures as much as FM 
majority spin. 


Exchange (meV) PBEO HSE06 HISS 
J c 4.0 3.3 4.2 

J ab 3.1 2.2 4.2 


TABLE II: Exchange interactions calculated within different 
approximations to the exchange-correlation functional for the 
crystal structure reported in Ref. 0- 


T (K) 3.5 60 90 100 120 160 200 300 

J c (meV) 3.0 2.2 5.8 3.7 9.6 3.4 4.8 2.7 

J ab (meV) 7.2 7.3 5.9 5.6 12.1 6.8 4.9 7.5 
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TABLE III: Exchange interactions calculated using HSE06 
hybrid-exchange functional for the crystal structures reported 
in Ref, [bill . 


2. Exchange interactions in CaCrOs, 

The computed exchange interactions for the two sets of 
crystal structures in Rehjis| and Ref.[3] have been tab¬ 
ulated in Table ITIl and Table Hill respectively. Although 
these exchange interactions are computed within differ¬ 
ent functionals, they have qualitatively the same struc¬ 
ture. The spins are coupled antiferromagnetically along 
the c-axis, but frustrated in the a6-plane as the NN 
exchange interactions is anti-ferromagnetic against the 
anti-ferro-orbtial ordering. According to Goodenough- 
Kanamori-Anderson (GKA) rules, this type of orbital 
ordering should result in a ferromagnetic spin coupling. 
However, the anti-ferromagnetic NN exchange interac¬ 
tion in the a6-plane frustrates this orbital ordering in the 
afe-plane, which might induce a large magnetoelastic cou¬ 
pling 0 . This could further be linked to the interplay 
between spin-ordering and lattice vibrations. 


B. The DFT + U electronic structure of CaCrOs 


FIG. 3: (Colour on line.) DFT band gaps in eV as a function 
of temperature are shown. The structures are adopted from 

Ref.fij. 


We have computed the exchange interactions J a & and 
J c and the energy difference between the AFM-A and 
AFM-C states, according to the DFT + U (with U = 1) 
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total energies. As shown in FigU below the Neel tem¬ 
perature, the energies of the AFM-A are higher than that 
of the AFM-C state. The Neel temperature is a turning 
point for the energies difference between the AFM-A and 
AFM-C states. Correspondingly, this point is important 
for the magnitude of exchange interactions. 



FIG. 4: (Colour on line.) The energy differences between 
AFM-C and AFM-A and exchange interactions computed by 
using DFT + U method (U = 1) are shown. The structures 
are adopted from Ref. 0 - 


C. Discussions 

1. Magnetic states 

The computed PDOS are qualitatively consistent with 
the x-ray absorption spectroscopy (XAS) as shown in 
Fig.2a of Ref.[|]. The two main peaks in the computed 
PDOS corroborate the main peak and the shoulder in the 
O 2 p intensity of the soft x-ray measurements in Ref. [§}. 
The measured energy gap between O 2p-orbital and d- 
orbital peaks is ~ 5 eV, which is slightly larger than 
the predicted (~ 3.7 eV). The accurate prediction of the 
key parameter, Hubbard- U, is necessary to understand 
the nature of the electronic state of CaCrOs, which can 
be performed by using DFT + U method. To improve 
further the theoretical description, more advanced com¬ 
putational methods, such as GW (Hj, are needed. The 
inclusion of the spin-orbit coupling (SOC) is another im¬ 
portant aspect However, these are beyond the scope of 
this paper. The Neel temperature T)v has been measured 
from the 1 /x versus T curve, which is ~ 90 K in Ref. @,0. 
This is in a good agreement with the computed exchange 
interactions. In addition, the experimentally observed 
magnetic moments which were measured to be 3.7hr 
Ref.Q, which is much larger than the computed (2.0 hb), 
probably owing to the unsaturated orbital component. 
However, the observed 1.2 /zb reported in Ref.[7] seems an 
indication of the existence of itinerant magnetism. The 
inconsistency between the measured magnetic moment 


could be attributed to the interplay between spin and 
orbital degrees of freedoms as mentioned previously. 


2. Electron-hole liquid 

So far the metallicity in CaCr 03 has predominantly 
been observed in optical measurements including optical 
reflectibility and conductivity and photoemissions sa 
From the calculations presented here, the FM band gap 
decreases sharply around T = 90 K, meanwhile the in¬ 
directness of the band gap has to involve phonon exci¬ 
tations during optical measurements. A closely related 
concept to the optical conductivity and band gap renor¬ 
malization is electron-hole liquid (EHL) [29), 1301 • EHL is 
a phenomena in which electron-hole density is sufficiently 
large to effect the electronic structure, i.e., the band gap 
renormalization. EHL has been observed in the mate¬ 
rials with cubic structure and indirect band gaps, and 
plays a significant role in reducing band gap, leading to 
a transition from insulator to semimetal [3l|. Recently, 
it has been shown that the band gap in BaTiOs nano¬ 
particles can be narrowed partially by EHL[32|. Com¬ 
bining with the small gap in the FM configuration as 
shown in Fig[3l at the Neel temperature the magnetic 
fluctuation might play an important role to the forma¬ 
tion of EHL in CaCrC> 3 . Based on the calculations pre¬ 
sented here and the previous experiments on EHL, we 
can speculate that the abnormal electronic properties in 
CaCrC >3 could be interpreted as an EHL assisted mainly 
by magnetic fluctuations, orbital ordering, and small FM 
band gap. The quasi half-metallicity due to the small FM 
band gap is crucial for the understanding of the electronic 
properties of CaCrC> 3 . 


3. Surface state 

In addition to the usefulness of EHL in plasmon en¬ 
hancements [HJ, the dangling bonds on the surface state 
in CaCrOs could be another source of metallicity. On the 
other hand, presumably the surface state of CaCrOs is 
conducting, this would be useful for exciting surface plas- 
mons. In conjugation with the spins carried by Cr 4+ , the 
CaCrOs could be a promising candidate for manipulation 
of spins through surface plasmons [34} . 


IV. CONCLUSIONS AND OUTLOOK 

In summary, the electronic structure and magnetic 
properties of Cr 4+ -based perovskites CaCrOs have been 
calculated within DFT by using hybrid-exchange den¬ 
sity functional PBEO, HSE06, and HISS. The computed 
PDOS is qualitatively in agreement with the previous soft 
x-ray measurements. The compute exchange interactions 
are consistent with the previous magnetic measurements. 
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Moreover, the spin densities have shown an anti-ferro- 
orbital ordering in the afr-plane. The analysis of the ex¬ 
change interactions and orbital ordering shows that there 
exists a spin-orbital frustration in the afo-plane, which 
could induce abnormal metallic state. The abrupt jump 
in the band gap of the FM majority spin between 90 K 
and 100 K corroborates the jump in the optical measure¬ 
ments reported previously. The PDOS have suggested a 
strong hybridisation between 3 d and 2 p orbitals, which 
agrees with the previous experimental works. It might be 
worthwhile performing calculations for the Hubbard-17, 
taking into account SOC and using more advanced com¬ 


putational methods to compute the DOS in the future, 
which would inspire more discussions about this fascinat¬ 
ing material. 
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